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Abstract. With Monte Carlo methods, we investigate the universality class of 
the depinning transition in the two-dimensional Ising model with quenched random 
fields. Based on the short-time dynamic approach, we accurately determine the 
depinning transition field and both static and dynamic critical exponents. The critical 
exponents vary significantly with the form and strength of the random fields, but 
exhibit independence on the updating schemes of the Monte Carlo algorithm. From 
the roughness exponents CiCloc an d ( s , one may judge that the depinning transition 
of the random-field Ising model belongs to the new dynamic universality class with 
C Cioc 7^ Cs and Cioc 7*= 1- The crossover from the second-order phase transition to 
the first-order one is observed for the uniform distribution of the random fields, but it 
is not present for the Gaussian distribution. 



PACS numbers: 64.60.Ht, 68.35.Rh, 05.10.Ln 
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1. Introduction 

In recent years many activities have been devoted to the study of dynamic processes far 
from equilibrium. An example is the domain-wall dynamics in disordered media, which is 
important from both experimental and theoretical perspectives. A crucial feature of the 
driven domain interface in disordered media is the depinning phase transition at zero 
temperature. For the domain-wall motion in ultrathin magnetic films, the Edwards- 
Wilkinson equation with quenched disorder (QEW) is a typical theoretical approach 
[H [21 [3j SI El El El E] - It is usually believed that different variants of the model belong 
to a same universality class. 

The QEW equation is a phenomenological model, without detailed microscopic 
structures and interactions of the materials. The domain wall in a two-dimensional 
system is effectively described by a single- valued elastic string. An alternative, possibly 
more realistic approach, may be to construct lattice models at the microscopic level. 
The two-dimensional random-field Ising model with a driving field (DRFIM) is the 
simplest example [91 [TQl [TTJ [121 H3]- I n t ne DRFIM model, overhangs and islands may 
be created during the dynamic evolution, and the domain wall is not single-valued 
and one dimensional. The DRFIM model is closer to experiments [T4"l [T5l [T6] and 
does not suffer from the theoretical self-inconsistence as in the QEW equation [IT]. 
Very recently it has been demonstrated that the DRFIM model may not belong to the 
universality class of the QEW equation [12l[T3]. To distinguish these two models, one 
needs accurate measurements of the critical exponents. In this respect, the short-time 
dynamic approach has been proven to be efficient [51 [HI [191 [2Ql (2U [121 EE] • 

The purpose of this paper is to clarify the universality class of the DRFIM model. 
In Refs. [121 H3], Monte Carlo simulations are performed at a fixed strength of the 
random fields with a uniform distribution, i.e., A = 1.5. In this paper, we systematically 
investigate the possible dependence of the critical exponents on the strength and form 
of the random fields, and the updating schemes of the Monte Carlo algorithm, including 
the crossover from the second-order transition to the first-order one. On the other hand, 
we measure not only the global and local roughness exponents ( and (i oc , but also the 
spectral roughness exponent £ s , to identify the dynamic universality class of the domain 
interface. Following Ref. [22], there are four different universality classes in the interface 
growth, namely, 



In Ref. [23], however, it suggests that there is a new universality class of anomalous 

roughening, with ( ^ Q oc ^ ( s and Q oc ^ 1. We verify that the DRFIM model belongs 
to such a new dynamic universality class. In section 2, the model and scaling analysis 
are described. In section 3, numerical results are presented. Section 4 includes the 
conclusions. 




Family- Vicsek 



(1) 




super-rough 
faceted. 
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2. Model and scaling analysis 

2.1. Model 

The two-dimensional (2D) random-field Ising model with a driving field is defined by 
the Hamiltonian 

n= - j£ s i s j -'£(h i + H)s i . (2) 

<ij> i 

where Si = ±1 is an Ising spin at site % of a rectangle lattice L x x L y , < ij > denotes a 
summation over nearest neighbors and if is a homogeneous driving field. The random- 
field hi may be distributed in different forms. A typical example is a uniform distribution 
within an interval [-A, A], following Refs. 0, Q21 EE]. We take A from 0.8J to 2.3J 
and set the coupling constant J = 1. Another example is a Gaussian distribution 
with mean zero and standard deviation a. Different form and strength of the random 
fields may lead to different kinds of phase transitions [2U [25]. At large A or a, the 
depinning transition is of second-order for both the uniform and Gaussian distributions. 
As A decreases, the depinning transition crosses over to a first-order one for the uniform 
distribution which is bounded. However, such a crossover is not present for the Gaussian 
distribution which is unbounded. 

We first concentrate on the uniform distribution of the random fields and the 
random-single-spin-flip Monte Carlo algorithm, and then extend the calculations to the 
Gaussian distribution and to different updating schemes. Simulations are performed 
at zero temperature with lattice sizes L y = 256, 512, 1024, and 2048 up to t max = 
1280,2560 or 5120 for different random fields. Total samples for average are about 
20000. Main results are presented with L y = 1024 and simulations of different L y 
confirm that finite-size effects are already negligibly small. L x is taken to be sufficiently 
large so that the boundary is not reached. Statistical errors are estimated by dividing the 
total samples in a few subgroups. If the fluctuation in the time direction is comparable 
with or larger than the statistical error, it will be taken into account. 

The initial state is a semi- ordered state with a perfect domain wall in the y direction. 
The periodic boundary condition is used in y direction, while the spins at the boundary 
in x direction are fixed. To eliminate the pinning effect irrelevant for disorder, we rotate 
the lattice such that the initial domain wall orients in the (11) direction of the lattice 
[9]. After preparing the initial state, the random- single- spin- flip Monte Carlo algorithm 
is performed. We randomly select one spin, and flip it if the total energy decreases 
after flipping, and if at least one of the nearest neighbors is already flipped, i.e., only 
spins at the interface are allowed to flip [251 12S] . A Monte Carlo time step is defined as 
L x x L y single spin flips. As time evolves, the domain wall propagates and roughens. In 
Fig. [H the snapshots of the dynamic evolution of the spin configuration are displayed 
for different random fields with the uniform distribution. As A increases, overhangs 
first appear then islands. 

Due to the existence of the overhangs and islands, there may be different ways 
to define the domain interface. Here, we adopt the definition with the magnetization. 
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Denoting a spin at site (x, y) by S xy (t), we introduce a line magnetization and height 
function of the domain interface [12l 



Lr=l 



(3) 



h(y,t) = ^-[m(y,t) + l}. (4) 
We then calculate the average velocity and roughness function of the domain interface, 

V (t) = d <y> , (5) 

W 2 (t) = (%,t) 2 )-(%,t)> 2 . (6) 

Here (• • •) represents both the statistical average and the average in y direction. 
A more informative quantity is the height correlation function, 

C(r,t) = ([h(y + r J t)-h(y,t)} 2 ). (7) 

It describes both the spatial correlation of the height function in y direction and the 
growth of the domain interface in x direction. Further, we consider the Fourier transform 
of the height function [221 123], 

h(k, t) = -L ^ [%, t) - (h(y, t))} exp{iky). (8) 

JLy y=l 



The structure factor is then defined, 

S(k,t) = (h(k,t)h{-k,t)). (9) 

To obtain the dynamic exponent z independently, we introduce an observable 

F(t) = [M {2 \t) - M(t) 2 }/oo 2 (t). (10) 

Here M{t) is the global magnetization and M^ 2 \t) is its second moment. Fit) is nothing 
but the ratio of the planar susceptibility and line susceptibility. 



2.2. Scaling analysis 

For the uniform distribution with A > 1 and Gaussian distribution of the random 
fields, the depinning transition is a second-order phase transition. The order parameter, 
i.e., the average velocity v(t) should obey the dynamic scaling theory supported by 
the renormalization-group calculations j28j [TH1 US]. In the critical regime, there are 
two spatial length scales in the dynamic relaxation, i.e., the nonequilibrium spatial 
correlation and the finite lattice size L y , scaling arguments lead to a dynamic 
scaling form [22 EEl EE] , 

v(t, r, L y ) = b-^ u G(b~ z t, b l ' v r, b~ l L y ). (11) 



Universality class of the depinning transition 



5 



Here, b is a rescaling factor, (3 and v are the static exponents, z is the dynamic exponent, 
and t = (H — H c )/H c . Taking b ~ t 1 /*, the dynamic scaling form is rewritten as 

v(t, r, L y ) = t-^ z G(l, t^r, t~ l / z L y ). (12) 

In the short-time regime, i.e., the regime with £(t) ~ t l / z <C L y , the finite-size effect is 
negligibly small, 

v(t,r) = t~^ uz G(t 1/uz r). (13) 

At the depinning transition point r = 0, a power-law behavior is obtained, 

v(t) ~ t- p/v *. (14) 

With Eq. (I13p . one may locate the transition field if c by searching for the best power-law 
behavior of v(t, r) [TBI EE]- To determine z/, one simply derives from Eq. ( ITBl) 

d r ln^,T)| T=0 ~i 1/l/ *. (15) 

In general, the roughness function w 2 (t) and height correlation function C(r,t) may not 
obey a perfect power-law behavior in early times. Due to the random updating scheme 
in numerical simulations, the domain interface and its velocity also roughen even without 
disorder (A = 0). This leads to rather strong correlations to scaling. To capture the 
dynamic effects of the disorder, we introduce the pure roughness function 

Duj 2 (t) = co 2 (t) -co 2 b (t), (16) 

and height correlation function 

DC(r,t) = C(r,t)-C b (r,t), (17) 

where oo b (t) and C b (r, t) are the roughness function and height correlation function for 
A = 0. For a large lattice L y and at the transition point H = H c , Du 2 (t) and DC(r,t) 
should obey the standard power-law scaling behaviors [29], [30] EH [12] , 

Doo 2 (t) ~ t 2 ^, (18) 

and 

V ' ; \ t 2 ^ z if < < r V ; 

Here £(£) ~ £ is the global roughness exponent, and (i oc is the local one. 
The structure factor should follow a scaling form [22], [23] 

S(k,t) = k-<X +r >f.(kt 1 '*), (20) 
where the scaling function takes the form 

f r \ S u2{( ' Cs) if n>>1 fort 
f'M~\u«» if „<1 ' (21) 

and £ s is the spectral roughness exponent. 

Since w 2 (£) describes the fluctuation in x direction and M^ 2 \t) — M{t) 2 includes 
those in both x and y directions, the dynamic exponent z can be determined 
independently by 

F(t) ~ m/Ly ~ t V 7^- (22) 
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3. Monte Carlo Simulation 

We present the results with the random-single-spin-flip Monte Carlo algorithm in the 
first two subsections, and those with different updating schemes in the third subsection. 

3.1. Uniform distribution of random fields 

In Fig. [2(a), the interface velocity v(t) is displayed for different driving fields and different 
strengths of the random fields. For example, at A = 2.0, it drops rapidly for smaller 
H, while approaches a constant for larger H. By searching for the best power-law 
behavior, one locates the transition fields H c = 1.2028(2) and 1.4599(2) for A = 1.3 
and 2.0, respectively. As shown in the figure, v(t) for both A = 1.3 and 2.0 at H c 
show almost perfect power-law behaviors starting from rather early times. According 
to Eq. (fi~4l) . one measures the exponent /3/vz = 0.244(2) and 0.199(2) from the slopes 
of the curves for A = 1.3 and 2.0, respectively. 

The transition field H c varies significantly with the strength of the random fields, 
as shown in Fig. [2(b), and all the measurements of H c are summarized in Table [U 
The results are qualitatively in agreement with the previous work [23], but with much 
better accuracy. For A > 1, the transition field H c < A, and the depinning transition 
is of second-order. The order parameter v(t) exhibits a nice power-law behavior, as 
shown in the Fig. [2(a). For A < 1, the velocity approaches a constant for H > A, 
while drops rapidly to zero once H is slightly smaller than A. This is a signal that 
the depinning transition is of first-order for A < 1, and the transition field H c = A. 
Following Fig. [2(b), the transition field H c increases with A, and can be fitted to an 
exponential function, i.e. H c = 1.70 — 2.04exp(— 1.08A). 

In Fig. [3(a), the dynamic observable F(t) in Eq. ( flOl) is plotted at H c . The 
power-law behavior is detected. According to Eq. (j22|) . 1/z — 0.718(7) and 0.788(5) 
are extracted respectively. To calculate the logarithmic derivative d T lnv(t,T) = 
d T v(t,r)/v(t,r), we interpolate v(t,r) between different H, e.g., in the interval 
[1.195,1.210] for A = 1.3, and [1.456,1.464] for A = 2.0. In Fig. [3(b), d T lnv(t,r) 
is plotted at H c . A power-law behavior is observed. Based on Eq. ( fT5{) . \jvz = 0.704(6) 
and 0.674(6) are derived from the slopes of the curves, respectively. 

In Fig. |4(a), the pure roughness function Du 2 (t) = u) 2 {t) — ulif) is plotted for 
A = 1.3 and 2.0 at the transition field H c . Here u) 2 (t) is the roughness function 
for A = 0. The curves for both A = 1.3 and 2.0 show almost perfect power-law 
behaviors from rather early times. According to Eq. ( fTB"l) . we measure the exponent 
2(/z = 1.676(6) and 1.665(9) from the slopes of the curves respectively. 

The structure factor S(k,t) is displayed for A = 2.0 at H c = 1.4599 in Fig. S(b). 
According to Eq. fl2Tl) . data collapse for different t is observed in the inset. The 
exponents used for the data collapse are ( = 1.06 and z = 1.27 obtained from Fig. [3(a) 
and 13(a) respectively. Obviously, for a small u = kt 1 ^ 1, u~( 2< * +1 ^ f s (u) approaches a 
constant. For a large u>l, u~^ +1 ^ f s (u) exhibits a power-law behavior. One extracts 
the exponent 2( s + 1 = 2.68(2) from the slope of the curves, i.e., ( s = 0.84(1). 
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In Fig. E(a), the pure height correlation function DC(r,t) is displayed for A = 2.0 
at H c = 1.4599. For a large r ^> e.g., r = 512, one extracts the exponent 

2(/z = 1.663(4) by Eq. ([I9~]) . consistent with that from Fig. [4(a) within errors. For a 
small r <C DC(r, t) should be independent of t for a normal interface with ( = Q oc , 
according to Eq. (fT9l) . But DC(r,t) of r = 2 clearly increases with time t. A power-law 
behavior is observed from rather early times. According to Eq. f fT9|) . the slope of the 
curve gives 2(£ — Cioc)/z — 0.657(4). From the measurements of ( and z from Fig. [3(a) 
and [4(a), we calculate the local roughness exponent (i oc = 0.643(6). 

The pure height correlation function DC(r,t) is plotted as a function of r in 
Fig. [5(b). For a small r <C £(£), there exists a power-law behavior, but this region is 
rather limited. For a large r, there emerge strong corrections to scaling. We introduce 
the form of corrections [29] , 

DC(r) ~ [tanh(r/c)] 2C '° c . (23) 

Here, c is an adjustable parameter. The fitting to the curve extends to a larger interval 
and yields (i oc = 0.645(6). It is in agreement with that from Fig. [5(a) within errors. 

From the measurements of /3/vz, 1/z, 1/vz, 2(/z, 2( s + 1, 2(( — Q oc )/z and 2Q oc 
in Fig. [2] to [5j respectively. We calculate the dynamic critical exponent z, the static 
critical exponents (3 and u, the roughness exponents (, (i oc and ( s . All the measurements 
are summarized in Table [fl compared with those of the QEW equation. In Fig. [6](a), 
the critical exponents z, is, and ( are displayed for different strengths of the random 
fields. The dynamic critical exponent z decreases with A, and the difference between 
the DRFIM model and QEW equation increases with A, up to 20 percent. The critical 
exponent v depends non-monotonously on A, and the difference between the DRFIM 
model and QEW equation reaches maximally for A = 1.3 and 1.5, over 20 percent. In 
Fig. [6(b), the critical exponent /3 decreases monotonously with A and reaches constant 
for A > 1.5. 

In Fig. 0(a), the roughness exponents (, (i oc and ( s are plotted for different strengths 
of the random fields. According to Ref. [22], the exponents C> 0oc an d C are n °t 
independent, and there are four different subclasses, namely, Family- Vicsek, super- 
rough, intrinsic and faceted. Following Fig. [7(a) and Tabled], ( = Q oc = ( s for A < 1, 
and it belongs to the Family- Vicsek universality class. In the crossover regime from 
the first-order transition to the second-order one, e.g., for A = 1.05 and 1.1, ( = ( s , 
the DRFIM model looks somewhat similar to the QEW equation, but there exists still 
a major difference, i.e., Q oc < 1 for the DRFIM model, while Q oc ~ 1 for the QEW 
equation. For A > 1.2, ( ^ Q oc ^ ( s and (i oc ^ 1, and it indicates that the DRFIM 
model belongs to the new universality class suggested in Ref. [23]. Additionally, one 
interestingly observes that the exponent ( s is in agreement with ((+(i oc )/2 within errors. 
The possible mechanism remains to be understood. 

For A > 1, the local roughness exponent Q oc decreases with A, and the difference 
between the DRFIM model and QEW equation increases with A, up to over 35 percent. 
In our opinions, it is mainly the overhangs and islands that induce this difference. 
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3.2. Gaussian distribution of random fields 

As shown in Fig. Efb), the transition field H c varies linearly with a for the Gaussian 
distribution of the random fields, and all the measurements of H c are summarized in 
Table |2j The results are in agreement with those in Ref. [21], but with higher accuracy. 
Here the depinning transition is always of second-order. The crossover from the second- 
order phase transition to the first-order one is not present, because the Gaussian 
distribution is unbounded. All the measurements of the critical exponents are also 
summarized in Table [21 In Fig. [6(a), variation of the critical exponents z, u, and £ with 
a for the Gaussian distribution of the random fields is displayed. The dynamic critical 
exponent z decreases monotonously with a, similar to that of the uniform distribution of 
the random fields. However, the static critical exponent v increases monotonously with 
a, different from that of the uniform distribution of the random fields. In our opinions, 
the non-monotonous behavior of the critical exponents for the uniform distribution of 
the random fields is induced by the crossover from the second-order transition to the 
first-order one. In Fig. [7(b), the monotonous decrease of the roughness exponents £, (i oc 
and ( s with a is displayed for the Gaussian distribution of the random fields. The results 
indicate ( ^ (i oc ^ ( s and Q oc ^ 1, the same as that of the uniform distribution of the 
random fields. 

As shown in Figs. [6] and [7J all the critical exponents z, v, (3, (, (i oc and ( s vary 
in principle with the strength of the random fields, and the dependence of the critical 
exponents on the strength of the random fields is also different for the uniform and 
Gaussian distributions. At the limit of the strong random fields, the critical exponents 
may tend to the same for the uniform and Gaussian distributions, but it is less clear for 
the roughness exponents Q oc and ( s . In other words, the strong universality is violated 
in the depinning transition. Both the crossover from the second-order transition to the 
first-order one in the uniform distribution of the random fields, and the dynamic effect 
of the overhangs and islands lead to this violation of the universality. It is pointed out in 
Ref. [13] that overhangs and islands play important roles in the domain-wall motion. As 
the strength of the random fields changes, the overhangs and islands alter their critical 
behavior and modify the critical exponents at the depinning transition. 

Although the strong universality is violated, we observe that the static exponent 
/3 looks somewhat 'universal'. For the Gaussian distribution of the random fields, /3 is 
almost independent of the strength of the random fields. For the uniform distribution 
of the random fields, (3 also reaches a constant value for A > 1.5 which coincides with 
that of the Gaussian distribution. Are there other combinations of the critical exponents 
which show a 'universal' behavior? In Fig. E](b) , variation of the exponents (2£+l) / z and 
(3/vz with the strength of the random fields is displayed. In particular, the exponent 
(2£ + l)/z looks 'universal', similar to the exponent /3. According to Eqs. ( I18jl and 
f )22|) . M^ 2 \t) — M(t) 2 ~ tC^+i)/ 2 . In other words, the scaling exponent of the planar 
susceptibility is 'universal'. Since M(t) is equivalent to the average height function, 
M^ 2 \t) — M(t) 2 is nothing but the fluctuation of the domain wall. 
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3.3. Updating schemes 

To further understand the universality at the depinning transition, we examine different 
update schemes in the Monte Carlo simulations, e.g., the sequential sweep and parallel 
sweep. For the sequential sweep, we select a spin sequentially by the row perpendicular 
to the domain wall, and flip it if the total energy decreases after flipping. For the 
parallel sweep, we divide the lattice into two 'checkerboard' sublattices. The spins 
in each sublattice do not directly interact each other. Therefore we may update all 
the spins in each sublattice simultaneously. Different updating schemes in the Monte 
Carlo algorithm do not modify the transition field and critical exponents, and all the 
measurements are summarized in Table [3j The small deviation in the transition field 
H c and the exponent (3 should be induced by corrections to scaling in the early times. 

Because of the presence of the random fields, the magnetic system has many 
local minima separated by sizeable energy barriers. During the dynamic evolution, 
the system may be trapped in a local minimum without finding the global one, as a 
consequence, remains far from equilibrium. The lifetime of a metastable state may be 
considered as infinite at zero-temperature. This type of behavior is well modeled by 
the single-spin-flip Monte Carlo dynamics [27]. In ordered to verify the reliability of 
the single-spin-flip Monte Carlo dynamics, the two-spin-flip algorithm has been also 
performed. The power-law behaviors in Fig [2] to |5] are observed. All the measurements 
of the depinning transition field and critical exponents are summarized in Table |3j The 
depinning transition field H c is changed, but all the exponents are in agreement within 
errors with those of the single-spin-flip dynamics [27] . 

4. Conclusion 

With Monte Carlo methods, we have simulated the dynamic relaxation of a domain 
wall in the two-dimensional DRFIM model. Based on the short-time dynamic scaling 
forms, the transition field and all critical exponents at the depinning transition are 
accurately determined. Since the measurements are carried out at the beginning of the 
time evolution, the dynamic approach does not suffer from critical slowing down. 

The results are first presented for the uniform distribution of the random fields with 
the single-spin-flip Monte Carlo dynamics. 

i) For A < 1, the depinning transition is of first-order, and the transition field 
H c = A. The exponents (3 and v do not exist, but the growth dynamics of the domain 
wall is still meaningful. Since ( = (i oc = ( s , it belongs to the Family- Vicsek universality 
class. 

ii) For A > 1, the depinning transition is of second-order, and the transition field 
H c < A. As shown in Fig. [6](a) and Fig. UJa), critical exponents vary significantly 
with A. Except for the dynamic exponent z, all other critical exponents indicate the 
crossover from the first-order transition to the second-order one. The exponent v and 
roughness exponents exhibit non-monotonous dependence on A. 
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iii) For A > 1 and A -> 1, e.g., at A = 1.05 and 1.1, the DRFIM model looks 
somewhat similar to the QEW equation, but Q oc ^ 1. If (i oc = 1, it would be in the so- 
called super-rough universality class as the QEW equation is. For A > 1.2, one observes 
clearly ( ^ Q oc ^ ( s , and Q oc < 1. It indicates that the DRFIM model belongs to the 
new dynamic universality class suggested in Ref. [23] . 

For comparison, simulations have been also performed for the Gaussian distribution 
of the random fields. The crossover from the second-order transition to the first-order 
one is not present in this case. Generally speaking, the strong universality is violated in 
the depinning transition, and the critical exponents vary significantly with the strength 
and form of the random fields. In our opinion, the difference between the uniform and 
Gaussian distributions of the random fields is mainly induced by the crossover from the 
second-order transition to the first-order one in the former case, and the dependence of 
the critical exponents on the strength of the random fields comes also from the overhangs 
and islands. Finally, our simulations confirm that the critical exponents are independent 
of the updating schemes. 
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v(t) 


F(t),C(r,t) 


u 2 (t),C(r,t) 


C(r,t) 


S(k,t) 


A 


H c 





V 


z 


C 


Cloc 


Cs 


< 1.0 


A 






1.50(1) 


0.49(1) 


0.48(2) 


0.49(1) 


1.05 


1.0447(1) 


0.621(8) 


1.35(3) 


1.49(1) 


1.00(1) 


0.86(1) 


1.02(2) 


1.1 


1.0833(2) 


0.411(4) 


1.12(2) 


1.46(1) 


1.10(1) 


0.83(1) 


1.10(1) 


1.2 


1.1485(2) 


0.388(3) 


1.08(2) 


1.44(1) 


1.19(1) 


0.80(1) 


1.01(1) 


1.3 


1.2028(2) 


0.347(3) 


1.02(2) 


1.39(1) 


1.17(1) 


0.78(1) 


0.99(1) 


1.5 


1.2933(2) 


0.295(3) 


1.02(3) 


1.33(1) 


1.14(1) 


0.74(1) 


0.94(1) 


1.7 


1.3670(2) 


0.296(3) 


1.10(2) 


1.30(1) 


1.10(1) 


0.69(1) 


0.91(1) 


2.0 


1.4599(2) 


0.295(3) 


1.17(2) 


1.27(1) 


1.06(1) 


0.65(1) 


0.84(1) 


2.3 


1.5398(2) 


0.299(3) 


1.21(2) 


1.26(1) 


1.04(1) 


0.61(1) 


0.82(1) 


QEW 




0.33(2) 


1.33(4) 


1.50(3) 


1.25(1) 


0.98(6) 


1.25(1) 



Table 1. The depinning transition field and critical exponents obtained for different 
strengths of the random fields with the uniform distribution are compared with those 
of the QEW equation in Refs. [H[5j|8]. 





v(t) 


F(t),C(r,t) 


io 2 (t),C(r,t) 


C(r,t) 


S(k,t) 


a 


H c 





V 


z 


C 


Cloc 


Cs 


0.8 


1.3607(3) 


0.334(3) 


1.07(2) 


1.36(1) 


1.15(1) 


0.75(1) 


1.03(1) 


1.0 


1.4205(3) 


0.306(3) 


1.09(2) 


1.34(1) 


1.14(1) 


0.73(1) 


0.94(1) 


1.2 


1.4779(3) 


0.293(3) 


1.11(2) 


1.31(1) 


1.10(1) 


0.68(1) 


0.91(1) 


1.4 


1.5355(3) 


0.291(3) 


1.15(2) 


1.29(1) 


1.07(1) 


0.64(1) 


0.89(1) 


1.8 


1.6505(3) 


0.296(3) 


1.17(2) 


1.27(1) 


1.04(1) 


0.59(1) 


0.78(1) 


2.2 


1.7647(3) 


0.295(3) 


1.19(2) 


1.24(1) 


1.00(1) 


0.54(1) 


0.70(1) 



Table 2. The depinning transition field and critical exponents for the Gaussian 
distribution of the random fields. 
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sweep 


H c v 


z 


c 


Cloc 


c. 


SSF 


random 
sequential 
parallel 


1.2028(2) 0.347(3) 1.02(2) 
1.2034(4) 0.345(6) 1.02(2) 
1.2033(4) 0.338(6) 1.00(2) 


1.39(1) 
1.40(1) 
1.40(1) 


1.17(1) 
1.17(1) 
1.17(1) 


0.78(1) 
0.77(1) 
0.77(1) 


0.99(1) 
0.99(1) 
0.99(1) 


TSF 


random 


0.7545(3) 0.340(6) 1.01(2) 


1.40(2) 


1.19(2) 


0.78(2) 


1.01(2) 



Table 3. The depinning transition field and critical exponents for the uniform 
distribution of the random fields at A = 1.3 with different update schemes. SSF: 
the single-spin-flip; TSF: the two-spin-flip. 




Figure 1. Snapshots of the spin configuration for the uniform distribution of the 
random fields. The black denotes the spin Si = 1 and the white denotes the spin 
Si = — 1. Simulations arc performed with the random-single-spin- flip Monte Carlo 
algorithm. 
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Figure 2. (a) Interface velocity v(t) is displayed for different driving fields H at 
A = 2.0 with solid lines for the uniform distribution of the random fields. For 
comparison, v(t) is also plotted for A = 1.3 and 2.0 at H c = 1.2028 and 1.4599 
respectively with circles and squares. Dashed lines show power-law fits. For clarity, 
the solid lines are shifted up. (b) The transition field H c is displayed, with A and a for 
the uniform and Gaussian distributions of the random fields respectively. The dotted 
and dashed lines show exponential and linear fits respectively. Errors are smaller than 
the symbol sizes. 
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Figure 3. (a) F(t) is plotted for A = 1.3 and 2.0 at H c . For clarity, the curve for 
A = 1.3 is shifted up. (b) The logarithmic derivative of v(t, r) is plotted for A = 1.3 
and 2.0 at H c . In both (a) and (b), the results are for the uniform distribution of the 
random fields, and dashed lines show power-law fits. 
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Figure 4. (a) The pure roughness function Dtu 2 (t) is displayed for A = 1.3 and 2.0 
at H c . (b) The structure factor S(k,t) is plotted for A = 2.0 at H c . In the inset, data 
collapse for different t is shown. In both (a) and (b), the results are for the uniform 
distribution of the random fields, and dashed lines show power-law fits. 
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Figure 5. (a) Time evolution of the pure height correlation function DC(r,t) is 
displayed for A = 2.0 at H c . Dashed lines show power-law fits, (b) DC(r, t) is plotted 
as a function of r for A = 2.0 at H c . The solid line shows a power-law fit with a 
hyperbolic tangent correction. In both (a) and (b), the results are for the uniform 
distribution of the random fields. 
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Figure 6. In (a) and (b), variation of the critical exponents z, v, £, (3, (2( + l)/z 
and P/uz with A and a respectively for the uniform and Gaussian distributions of the 
random fields is displayed. 
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Figure 7. (a) The roughness exponents £, Ooc, Cs and (C + 0oc)/2 are plotted for the 
uniform distribution of the random hclds. (b) The roughness exponents £, Cloct Cs and 
(C + 0oc)/2 are displayed for the Gaussian distribution of the random fields. In both 
(a) and (b), full lines are just to guide the eyes. 



